UTTG-04-08 



A Tree Theorem for Inflation 



Steven Weinberg* 
Theory Group, Department of Physics, University of Texas 
Austin, TX, 18112 



Abstract 
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I. Introduction 



Quantum field theory is used differently in cosmology than in elementary 
particle physics. In particle physics we need to know matrix elements be- 
tween "in" states and "out" states, defined respectively by their appearance 
at times long before and long after a collision. In contrast, in cosmology we 
are interested in the expectation value of various operators in an "in" state, 
usually defined as a state that looks like the vacuum at very early times, 
long before perturbations leave the horizon during inflation. To calculate 
such expectation values, we must use the "in-in" formalism of Schwinger et 
al.[l] 

The "in-in" formalism provides graphical rules for calculating expecta- 
tion values, similar to but more complicated than the Feynman rules used 
in particle physics. A special role in these calculations is provided by tree 
graphs, which in many theories give much larger contributions to expecta- 
tion values than graphs with loops. For instance, in general relativity the 
coupling constant G appears only in a factor 1/G multiplying the whole 
Einstein-Hilbert action, so each graviton propagator yields a factor G and 
each interaction comes with a factor 1 JG. A connected graph with / internal 
lines and V vertices thus yields a factor G J ~ V = G L_1 , where L = I — V +1 
is the number of loops, so at a characteristic frequency H, graphs with L 
loops are suppressed relative to tree graphs by factors of order (GH 2 ) L . The 
strength of fluctuations in the cosmic microwave background suggests that 
GH 2 ~ 1CP 10 . This suppression can partly be compensated by the appear- 
ance of powers of In a [2] (where a is the Robertson- Walker scale factor), but 
for this to make loops competitive with trees the universe would have to 
expand after horizon exit by something like 10 10 e-foldings. For the same 
reason, tree graphs will dominate in any theory with a small very coupling 
constant g that appears only as a factor 1/g in the action. 

In the present state of cosmology it is fortunate that tree graphs make 
a much larger contribution to correlation functions than graphs with loops. 
During the long period from when perturbations left the horizon during 
inflation to when they re-entered the horizon in the radiation- or matter- 
dominated era, there were various events, such as reheating, lepton and 
baryon synthesis, and dark matter decoupling, about which we know es- 
sentially nothing. The only reason that we are able to relate observations 
of the microwave background anisotropics or large scale structure to what 
happened in inflation is that certain quantities such as the curvature per- 
turbation ( and the gravitational wave amplitude are believed to be time- 
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independent when perturbations are outside the horizon[3] — that is, when 
the physical wave number k/a\s much less than the expansion rate H. For 
tree graphs, when the physical wave numbers associated with external lines 
are all much less than H then the same is true of the wave numbers asso- 
ciated with all internal lines but, as shown by the In a factors mentioned 
above, this is not true of graphs with loops, and so the theorems that state 
the constancy of Q or of the gravitational wave amplitude outside the horizon 
apply only to tree graphs. 

It is well known in elementary particle physics that the tree graph con- 
tributions to matrix elements between "in" and "out" states can be calcu- 
lated by solving the classical field equations in the presence of an external 
c-number current. [4] In the Appendix to a recent paper [5] I remarked that 
the same is true in the "in-in" formalism, but the prescription stated there 
was not correct. Section II of the present paper states and proves a general 
theorem (hopefully correct this time) for the evaluation of tree contributions 
to expectation values by the use of classical field equations. In Sections III 
and IV this theorem is applied to inflation. The derivation of the tree the- 
orem in Section II relies on an analytic continuation of trajectories in the 
path-integral formalism, and since this is not treated rigorously, we check in 
Section V that this theorem does give the correct results in the first few or- 
ders of perturbation theory for the inflationary model discussed in Sections 
III and IV. The methods of this paper are not really needed in the calculation 
of tree contributions to specific expectation values, since the rules for such 
calculations are already well known, but it is hoped that our results may 
prove useful in proving general theorems about non-Gaussian correlations 
in cosmology. 

II. Tree Theorem for the "In-In" Formalism 

We will consider a general Hamiltonian system, with Hermitian operators 
Q a (t) and their canonical conjugates P a (t), satisfying the usual commutation 
relations 

[Q a (t),P b (t)} = i5 ab , [Q a (t), Q b (t)} = [P a (t), P b (t)} = . (1) 

In field theories a is a compound index, including a spatial coordinate as 
well as discrete indices labeling the nature of the various operators, and 
5 ab includes a delta function in the spatial coordinates as well as Kronecker 
deltas for the discrete indices. These are Heisenberg-picture operators, with 
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a time-dependence generated by a Hamiltonian H[Q(t), P(t),t]: 



Qa(t) = -i[Q a (t),H[Q(t),P(t),t]\ , P a (t) = -i[P a (t),H[Q(t),P(t),t] 

(2) 

We are allowing the Hamiltonian here to have an explicit dependence on time 
for a reason discussed in [5]: When Q a (t) and P a (t) are the fluctuations of 
canonical variables around time-dependent background values, the Hamil- 
tonian that generates their time-dependence is not the time-independent 
Hamiltonian for the total canonical variables, but the sum of the terms in 
this Hamiltonian of second and higher order in the fluctuations, so that a 
time-dependence is introduced by dropping the terms of zeroth and first or- 
der in fluctuations. The differential equations (2) show that in H[Q(t), P(t),t] 
we can set the time argument of Q a (t) and P a (t) equal to any fixed time i*, 
but the explicit time dependence of H[Q(t#), P(t*),t] still remains. 

Instead of giving a formula for expectation values of products of the 
Qa{ti) at a fixed time t\ in an "in" state |0,in) defined by its appearance 
at a time to < ti, we will give a formula for a generating function W[J, ti] 
of a c-number current J a , from which all such expectation values may be 
obtained. The generating function is defined by 



e w[JM= o , in 



exp ^2Q a {h)J a 



0,in 



(3) 



and from its derivatives we can calculate the expectation value of a product 
of n of the Q's: 



<0,in|Q a (ti)Q 6 (ti)--.|0,in) = 



dJ a dJ b 



exp 



[W[j,h] 



(4) 



(In field theories, sums over a like that in (3) include an integral over a 
spatial coordinate as well as sums over discrete indices, and the derivatives 
in (4) are functional derivatives.) We need not take J a to be real, but if we 
do then the generating function is real also. 

To calculate W[J, t±] in the tree approximation we introduce a pair of 
complex c-number J-dependent functions of time, qLa{t) and qn a {t), which 
are defined by three conditions: 



(A) Both qha{t) and qn a {t) satisfy the Lagrangian equations of motion 

. dL[q L (t),q L {t),t] 



d_ / dL[q L (t),q L (t),t] 
dt V dq La (t) 



dqLa{t) 



(5) 



4 



and 

d fdL[q R (t),q R (t),t] \ _ dL[q R (t),q R (t),t] 



dt V dq Ra (t) ) dq Ra (t) (6) 

where L[q(t),q(t),t] is obtained from the classical Hamiltonian H by using 
the expression 

L[q(t),q(t),t] = J2 Pa (t)q a (t) - H[q(t),p(t),t] . (7) 

a 

with p a {t) eliminated from the right-hand side by using the classical formula 

_ dH[q(t),p(t),t] 
Qait) ~ d Pa (t) ■ (8) 

(B) The qLa{t) and q Ra (t) and their time derivatives satisfy constraints at 
the time t±: 

QLaih) = QRaih) , (9) 

and 

dL[q L (ti), qiih), h] dL[q R {t{), q R (h), h] 



dqLa(h) dq Ra (ti) 



-iJa ■ (10) 



(C) The qLa{t) and qRa(t) and their time derivatives also satisfy constraints 
at the time to that is used to define the state |0, in), constraints that depend 
on the nature of this state. In particular, if |0, in) is a state that looks like 
the Bunch-Davies vacuum at a time to = — oo, then qLa(t) and q Ra (t) satisfy 
"positive frequency" and "negative frequency" conditions, respectively; that 
is, they must for t — > — oo be superpositions of terms with time-dependence 
proportional to eT %u)t and e +lujt , respectively, where the cjs are various posi- 
tive frequencies. 

With qha{t) and q Ra {t) calculated in terms of the current using these 
three conditions, the generating function W[J, t\] is given in the tree ap- 
proximation by 

W[J, ii]tree = i t dt {L[q R (t),q R (t),t] - L[q L (t), q L (t), t]\ +]T q Ra (ti) J a ■ 

J to a 

(11) 

(There is an easy generalization of this theorem: Instead of the linear 
function J2 a Qa(ti)J a hi the exponential on the right-hand side of the defi- 
nition (3) of W, we could insert an arbitrary function J[Q{ti)\ of all the Q a 
at the same time t\. Then in the tree approximation W would be given by 

W tvcc = ij^ dt {L[q R (t),q R (t),t} - L[q L (t),q L (t),t}} + J[q R (h)} 
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with J a on the right-hand side of the constraint (10) replaced with dj(q)/dq a , 
in which we set q a = qR a (h).) 

The proof of the tree theorem relies on the path-integral formulation 
of the "in-in" formalism, so we begin with a brief derivation of the path- 
integral formula for exp(W^[J, ti]). At any given time t, we introduce a 
complete set of states \q, t), defined as normalized eigenstates of the Q a (t) 
with eigenvalues q a : 

Qa(t)\q,t) = q a \q,t) , {q, t\q f , t) = J] 6(q a - q' a ) . (12) 

a 

The expectation value (3) may then be written 

e W[J ' tl] = J (j\dq a dg' a d^^)(q,t 1 W,t r 

X exp (^2 QaJa^j (q, h\q", t >*o(«") , (13) 

where ^o(q) is the wave function of the state in which the expectation value 
is taken 

*o(g) = <g,*o|0,in) . (14) 

The matrix elements between eigenstates of the Q's at times t\ and to with 
to < t\ are given by an integral over real functions q a (t) that interpolate be- 
tween the eigenvalues at to and t\ together with independent unconstrained 
real functions p a (t)- 

(q,h\q',t ) = J \J[dq a (t)dp a (t)) (jJS(q a (t ) - q'j'j (jl <%a(*l) - <Za)) 



x exp I i dt 



J2p a (t)q a (t) ~ H[q(t),p(t),t] 



(15) 



(Eq. (15) and its derivation are the same as encountered in the familiar 
derivation of the path- integral formula for the S-matrix.) Using this in 
Eq. (13), we must introduce separate integration variables qLa(t), PLa{t) and 
QRa(t), PRa(t) in the path-integral formulas for {q,ti\q',to} and {q, ti\q", to), 
respectively. This gives 

e W[J, tl ] = j mdq La (t)d PLa (t)) j \l[dq Ra (t)dp Ra (t) 
J \a,t I J \a,t 



X 



X 



*o(«l(^)*o(Wo)) [jlHqLaih) -qRa{h))J exp (^2q Ra (ti)Jaj 

exp [^ij^dt \^VLa(t)q L a(t) - H[q L {t),p L (t),t]^j 

exp (ij* 1 dt h2 PRa (t)q Ra (t) - H[q R (t),p R (t),t]V\ . (16) 



There is no special reason why we chose q Ra (ti) rather than qLa(t\) to multi- 
ply the current in the first exponential; the delta function makes this choice 
inconsequential. 

Eq. (16) leads to well-known graphical rules: Writing H as the sum of 
a quadratic part and an interaction, we expand in powers of the interaction 
and J, and evaluate the resulting Gaussian integral as a sum of ways of 
pairing the qs and ps in this expansion. Because it enters in the exponential 
on the left-hand side of Eq.(16), W[J, ti] is given by a sum of connected 
graphs. To isolate the connected tree graphs, we use the same trick as is 
used in proving a tree theorem for the S-matrix[4]: We introduce a fictitious 
coupling constant g, defining a generating function W[J, h,g] by 



e W ^]l9^ J \l[dq La (t)dp La (t)j J \l[dq Ra (t)dp Ra (t) 

X *o(^(*o))* (Wo)) (lI%La(il) " 9Jta(*l))) 
x ex P ^X^fla(*i) J aj exp (y~ l -J t dt \^ZPLa{t)qLa{t) - H [q L (t), p L (t), t] j^j 
x exp (^J^dt i^2pRa{t)QRa{t) - H[q R (t),p R (t),t]^j . (17) 

By the same argument as in Section I, a connected graph with L loops makes 
a contribution to W[J,t\,g]/g proportional to g L_1 , so W[J,ti,g] for g — > 
approaches a g-independent limit given by the sum of tree graphs for any 
value of g, and in particular for the physical value g = 1: 

W[J,t 1 ,g]^W[J,t 1 ] tree for g^O. (18) 

Now, in the limit g — > 0, the path integrals in Eq. (17) are dominated 
by complex J-dependent trajectories of qLa(t), pia{t) and q Ra (t), piu^t) at 
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which the coefficient of 1/g in the combined exponential is stationary, so 



W[J,t 



1 J tree 



E qRa{tl)Ja ~ i f 1 dt \Y,PLa{t)qLa(t) ~ H[q L (t) , p L (t) , t] 
+i dt | E PBa {t)qna (t) -#[<?«(*), p fi (t) , t] 



with the condition that the trajectories be stationary (indicated by the sub- 
script "staty") taken subject to the constraint (9) imposed by the factor 
Eta fi(lLa(ti) — QRa(ti)) in Eq. (17), and also subject to constraints imposed 
by the factors ^o^l^o)) an d (?i?(*o)) j about which more later. 

The functions PLa(t) and pR a (t) are unconstrained, so in Eq. (19) we 
may take them to be given by the classical conditions (8) applied to both 
functions. Then Eq. (19) becomes 



W[J,tl]tree 



rti 

22qRa(h)Ja-i dt L[q L (t),q L (t),t] 

a Jt 
rti 

+i dt L[q R (t),q R (t),t] 
Jt 



(20) 



with L[q,q,t] given by using (8) to eliminate the ps in Eq. (7). This is the 
same as the desired result (11), provided we can show that the trajectories 
for which this quantity is stationary are those described by the three above 
conditions (A), (B), and (C). 

To implement the condition that this is stationary with respect to vari- 
ations in qL a (t) an d QRa(t), we note that 

J2<lMtl)Ja-i I 1 dt L[q L (t),q L (t),t]+i f ' dt L[q R (t),q R (t),t] 

n J to Jtn 



Jtn „ " " ' 



'to a - 

rti 

+ i 

I to 



1*? 



dqia{t) 
dL[q R (t),q R {t),t] 



^Sqitaih) (,J a + 



dqRa{t) 
. dL[q R (ti),q R {ti),ti] 
dq R a(h) 



dqiait) 

dL[q R (t),q R (t),t] x . . 
5qRa{t) H tt. — -t 5q Ra {t) 



dqRa{t) 
'^SqLaih) 



dL[q L (ti),q L (ti),ti 
dqLaih) 



8 



: / dt^ 



dL[q L (t),q L (t),t] d dL[q L (t),q L (t),t] 



dq La {t) dt dq a {t) 



SqLa(t) 



+ i 



— WL[q R (t),q R (t),t] d dL[q R (t),q R (t),ty gi ir {r>]) 



J my. 

Jto a 



dq Ra (t) dt dq Ra (t) 



The vanishing of the coefficients of SqL a (t) and Sq Ra (t) for to < t < t\ 
yields the Lagrangian equations of motion (5) and (6) . The vanishing of the 
coefficients of 5qLa{ti) and Sq Ra (ti), subject to the condition that 5qLa{ti) = 
<W*i), yields Eq. (10). 

Finally, we must return to the constraint on the solutions of the La- 
grangian equations of motion (5) and (6) imposed by the presence in (17) 
of the wave functions ^o(^q R (to)^ an d ^(^(^o)) . As is familiar from the 
path-integral calculation of the S-matrix, where to ^ —oo and |0,in) is the 
"in" vacuum the factor ^o(q R (to)^J has the effect of putting a —it in the 
denominator of the Fourier integral for the propagator, so when Eq. (6) 
is solved using this propagator as a Green's function we get "negative fre- 
quency" solutions q Ra (t) — that is, solutions that behave for t — > — oo as 
a superposition of terms with time dependence e tujt , where to > 0. As re- 
marked in [5], because L[qL(ti),qL(ti),ti] enters in the argument of the 
exponential in (17) with an opposite sign to L[q R (ti),q R (t{),ti], the effect 
of the wave function ^(^(^o)) is to constrain qLa(t) to have "positive fre- 
quency" — that is, to behave for t — > — oo as a superposition of terms with 
time dependence e~' lujt , again with lo > 0. 

This completes the proof that the functions qLa{t) and q Ra (t) that should 
be used in Eq. (11) to calculate the tree contribution to W[J, t±] are those 
satisfying the above conditions (A), (B), and (C). In the case where J a is real, 
these three conditions are consistent with the result that qLa(t)* = qRa(t), 
which makes the tree approximation (11) to the generating function W[J, ti] 
real, as expected. 

III. Application to Inflation 

To illustrate the use of the tree theorem of Section II in cosmology, we will 
consider a simple semi-realistic model, in which a single scalar field evolves in 
an unperturbed Robertson- Walker metric g^ u of zero spatial curvature. We 
write the scalar field 0(x, t) as an unperturbed term (p(t) plus a fluctuation 
<p(x.,t) (called Sip(x, t) in [5]). The Lagrangian for the fluctuation is 



L = J Vgd 3 x 
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°7 



(22) 



where a(t) is the Robertson-Walker scale factor, and in accordance with the 
prescription for constructing the Hamiltonian for fluctuations mentioned in 
Section II, the time-dependent potential U (99, t) for the fluctuations is a 
function containing only second and higher powers of (p{x,t), and given in 
terms of the time-independent potential V(4>) for the total scalar field by 

U(<p,t) = v($(t) + <p) - Vfatj) - V'(4(tj)<p . (23) 

It is straightforward to apply the results of the previous section to this 
theory. Here the index a is the spatial coordinate x; the variable q a (t) is 
(p(x,t); the current J a is J(x); and derivatives with respect to q a (t) or q a (t) 
are functional derivatives. The Lagrangian equations (5) and (6) are here 
the Euler-Lagrange equations 

q> L + 3Hip L - a~ 2 V 2 ip L + U'(<p L ,t) = , (24) 

<PR + 3H<p R -a- 2 V 2 <p R + U'(<p R ,t) = 0, (25) 

where H(t) = a(t)/a(t), the prime on U indicates a derivative with respect 
to its field argument, and the constraints (9) and (10) read 

¥>l(x,<i) = ipR(-x.,h) , (26) 

L (x, h) - <p R (x, h) = -io _3 (ti) J(x) . (27) 

With to = —00 and expectation values calculated for the Bunch-Davies vac- 
uum, the "positive frequency" and "negative frequency" constraints require 
in this model that, for t — > —00, <^l(x, t) and </?r(x, t) approach superposi- 
tions of exp(ik • x — ikrf) and exp(ik • x + ikrj), respectively, where r\ is the 
conformal time, with 77 > 0. 

After solving Eqs. (24) and (25) subject to these constraints, the gener- 
ating functional W[J,ti] can be calculated in the tree approximation from 
Eq. (11), which for this model reads 

W[J,h] tree = i J' dt{L[(pR(t),(pR(t),t]-L[<p L (t),(p L {t),t]}+J dFxtpRfch) 

(28) 
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with the Lagrangian L given by Eq. (22). The integral of the Lagrangian 
can be simplified by integrating by parts and then using the field equation. 
Note that 



j ' 1 dt J d 3 x a 3 (t) 
- [ 1 dt f d 3 x a 3 (t) 

J — OC J 



1 



a 3 (tiMtiM£i 



1 



cp{ if + 3Hp - ^ VV) + U((p, t) 



(There is no contribution from the lower limit of the integral, because the 
integrand oscillates increasingly rapidly for t — > — oo.) Hence by using the 
field equations (24) and (25), Eq. (28) becomes 

W[J,tl]tree = ^ (ti) j d 3 X {<£fl(x, ti)(f> R (x., ti) - </?l(x, ti)(f> L (x., ti)} 

+i p dta 3 {t) j d 3 x f^(p R (x,t)U'{ip R {x,t),t) -u{ip R (x,t),t) 

- ^ L (x, t) 17' (^(x, t), t) + £/(Wx, i), t) | 

+ y a! 3 x ¥> fl (x,ti) J(x) . 

Using the constraints (26) and (27), we see that the first term is —1/2 the 
last term, so 

W[J,h] tree = i p dta 3 (t) J d 3 x f^<p R (x,t)U'(<p R (x,t),t) -u{ip R (x, 

- ^ L (x, t)C/'(Wx, i), t) + 17 (^(x, i), *) | 

+ i /" d 3 x ¥> fl (x,ti) J(x) . 

This is the form we will use in what follows. With W calculated in this 
way, the tree approximation to the expectation value of a product of n ips 
is given by Eq. (4) as 

6 n 



(0, in |y>(x,ti)<p(y,ti) ■■■!(), in) 



tree 



cxp 



<5J(x)<5J(y) 
IV. Integral Equation Formulation 



W[J,h] tiee 



j=o 
(30) 
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We will now re-write the Euler-Lagrange equations as integral equations 
that incorporate the constraints on the behavior of </?l(x, t) and </?r(x, t) for 
t — > — oo as well as the constraints (26) and (27). It is these integral equa- 
tions that will be used in the following section to check that this formalism 
generates the usual results of perturbation theory. 

We first separate the potential U into a term quadratic in ip and an 
interaction term T containing only cubic and higher terms: 

U(<p,t) = ±<p 2 U"(0,t) + r(<p,t) . (31) 

We introduce functions Uk{t) for which u k (t) exp(ik • x) are "positive fre- 
quency" solutions of the linearized Euler-Lagrange equations, — that is, 

u k {t) + 3H(t)u k (t) + (k/a{t)fu k (t) + U"(0,t)u k (t) = . (32) 

with "positive frequency" interpreted to mean that the WKB solution for 
t — > — oo is proportional to 



a(t)" 1 exp (^-ik J dt'/a{t')j , 



rather than to its complex conjugate. It will be convenient to normalize 
these functions so that for t — ► — oo the WKB solution takes the form 

Uk{t) ^ (2vr)3/ 2v ^ a( t) eXP {- iH fuW)) ' ^ 

with t* any fixed time. The Wronskian of u k and u* k is proportional for all 
times to 1/a 3 , so the normalization (33) gives it the value 

1 

U k (t)u* k (t) - U k (t)u* k (t) = ^3^) • ( 34 ) 

Using these functions, we can construct a Green's function: 

G(x-x',M') = 9(t-t')(G (x-x',t,t') + G* (x:-x.',t,t'f) . (35) 
where 6 is the usual step function, and 

G (x - x , t, t') = -if d 3 k exp (ik • (x - x')) u k {t)u%{t') . (36) 
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Both Go and its complex conjugate satisfy the homogeneous wave equation 
+ 3H ^]f t ~ «~ 2 (*)V 2 + U"(0, t)j G (x - x', t, t') = (37) 

and the Wronskian formula (34) then tells us that G satisfies the correspond- 
ing inhomogeneous equation: 

+ 3H(t)^ - a- 2 (t)V 2 + U"(0, t^j G(x-x', t, t') = -6(t-t')6 3 (x-x') . 

(38) 

Solutions of Eqs. (24) and (25) that satisfy the constraints (26) and (27) 
and the positive and negative frequency conditions, respectively, are given 
by the coupled integral equations 



x,t) = J d 3 x' p dt' S^G(x-x',t,t')a 3 (t')T'(ip L (x',t'),t^ 



+ G (x-x',t,t')a 3 (t^ 



+ i J d 3 x' G (x - x , t, h) J(x') , 



(39) 



<p R {x,t) = J^x'j' 1 dt' ^G(x-x',t,t')a 3 (t')r'(cp R (x',t'),t' 



G^x-x'^W 



r'^ R ( x ',t'),t')-r'^ L (x',t'),t' 



-i J d 3 x' GS(x-x',t,ti) J(x') . 



(40) 



Functions y>£,(x,t) and </?r(x, t) that satisfy these integral equations will 
satisfy the field equations (24) and (25) because of the properties (37) and 
(38) of the Green's functions, and they will satisfy the positive and negative 
frequency conditions because Go and Gq respectively satisfy these condi- 
tions, while G(x — x',i,i') — > for t — > — oo. To check the constraint (26), 
and for future reference, we note that since 9{t\ — if) = 1, 

99 L (x,ti) = <£>fl(x,ii) 

= j d z x' ^ / |G (x-x / ,t 1 ,t , )a 3 (t , )r , (^(x / ,t / ),i / ) 
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+ G^(x - x', t u t') a 3 (t') r'(^(x', t'),t') 
+ f dVD(x-x',ti) J(x') , 



(41) 



where D is the real function 

D(x — x', t\) = iGo(x — x', ti, ti) = — iGg(x — x', t\,t\) 

= J d 3 k \u k (h)\ 2 exp (ik • (x - x')) • (42) 

To check the constraint (27), we note that Go(x — x', t,t) is imaginary, so 

G(x - x', t, t') = 6(t - t') (g (x - x', t, t') + G* (x - x', t, 0) , (43) 

the dot indicating a derivative with respect to the first time argument. It 
follows then from Eqs. (39) and (40) that 

<f L (x,t 1 )-(p R (x,t 1 ) = i J dV (Go(x-x^i,*i) + Go(x-x',ii,ii)) J(x') . 

(44) 

Using the Wronskian formula (34) again gives 

G (x - x', h, h) + GS(x - x', ti, ti) = -a" 3 (t!) ,5 s (x - x') , (45) 

so Eq. (44) shows that the constraint (27) also is satisfied. 

For a real current J(x), and a potential U(<p) that satisfies the real- 
ity condition U*(<p) = U(<p*), the iterative solution of Eqs. (39) and (40) 
obviously satisfies the condition <£>l(x, t)* = <pr(x, t). For t = ti, when 
</?l(x, t±) = (Pr(x, ti), this implies that <^i?( x ) h) is real. It then follows that 
the formula (29) gives W[J, t±] real in the tree approximation. 

V. Perturbation Theory 

The derivation of the tree theorem in Section II was something less than 
mathematically rigorous. In particular, in evaluating the path integral in the 
limit g — > by setting the argument of the exponential equal to its value 
where stationary, we distorted the contour of integration of the variables 
q a (t) and p a (t) away from the real axis, without proving the validity of this 
analytic continuation. As a check on the validity of this theorem, in this 
section we will verify that when we expand the integral equations derived in 
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the previous section in powers of the interaction T, we obtain the first few 
orders of perturbation theory given by the usual path-integral formulation 
of the "in-in" formalism. 

Let us first recall the diagrammatic rules given by perturbation theory 
for the expectation value 



(0,in|p(x,ti)y>(y,*i)---|0,in) , 



(46) 



(as derived for instance in [5]), in the special case of the model discussed 
in Sections III and IV. We sum contributions from diagrams containing any 
number of vertices, each of which can be either of L or R type, and is labeled 
with a spatial coordinate and a time t<t\. For each vertex of L or R type 
we include a factor +ia 3 (t) or —ia 3 (t), respectively, together with a time- 
dependent coupling constant equal to the n-th derivative with respect to 
<p of F((p,i) at <f = for a vertex to which is attached n internal and/or 
external lines. Attached to the vertices are external lines, corresponding to 
fields in the expectation value (46), and/or internal lines connecting pairs 
of vertices. The contribution of an external line corresponding to a field 
</?(x, ti) attached to a L or R vertex labeled with spacetime coordinates y, 
t, is a function 

L: (0|/(y,t)/(x,t 1 )|0) (47) 



or 



R: ((^(/(^(/(y^lO) , 
respectively, where </? 7 (x, t) is the interaction picture field 

/(x,t)= ( d 3 k \e ik -*u k (t)a(k) +e- ik ' x 4(t)a*(k) 



(48) 



(49) 



with a(k) and a*(k) the usual annihilation and creation operators, and |0) 
is the bare vacuum, annihilated by a(k). The contribution of an internal 
line connecting two vertices of L and/or R type labeled with spacetime 
coordinates x, t and x', t' is 



LL 
RR 
LR 
RL 



(0|T{/(x,t)/(x',t')}|0), 
(0|T{/(x,t)/(x',t')}|0), 
(0|/(x,t)/(x',O|0), 
<(V(x',tV(x,t)|0), 



(50) 
(51) 
(52) 
(53) 



where T and T denote time-ordered and anti-time-ordered products, respec- 
tively. In addition, there may be lines passing through the diagram without 
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interaction, corresponding to pairings of fields in the product whose expec- 
tation value is being calculated. Each such line, corresponding to fields 
¥>(x,ii) and ip(y,t\), makes a contribution 

<(V(x,ti)/(y,ti)|0> , (54) 

the order of the fields here being unimportant, since fields at the same time 
commute. Finally, we are to integrate the spatial coordinates associated with 
each vertex over all space, and integrate the corresponding time coordinates 
from — oo to t±. 

As a first step in checking the agreement of our theorem with these 
graphical rules, we note that the quadratic term in U(<p) drops out in the 
combination ^ipU'(ip,t) — U(ip,t), so Eq. (29) may be written in terms of 
the interaction T defined by Eq. (31): 

w[j,h] trcc = i J' 1 dt J d 3 x |^i?(x,t)r'(^(x,t),t) - r(<p R (x,t),t) 

-\pL (x, t)V (ip L (x, t) , t) + r (ifL (x, t), t) | a 3 (t) 

+^Jd 3 x^ R (x,t 1 )J(x). (55) 

We can now expand in powers of V. 
Zeroth Order 

The zeroth order approximation to <^h(x, ti) is given by dropping the 
terms involving V in Eq. (41): 

(p<® ( x , ti) = J d 3 x L>(x - x', ti) J(x') . (56) 

Using this in Eq. (55) gives the zeroth order perturbation to the generating 
function: 

[J, h] tiee = ^Jd 3 xJ d 3 x> £>(x - x', h) J(x) J(x') . (57) 

The result that this is quadratic in J corresponds to the elementary obser- 
vation that in the absence of interactions all fluctuations are Gaussian. The 
two-point correlation is given in the tree approximation by Eqs. (30) and 
(57) as 

(0,in|^(x,ti)v9(y,ti)|0,in) trcc = D(x-y,ii) . (58) 
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This agrees with the above perturbative rules, because the two-point func- 
tion (54) is just D(x — y, t±). 

First Order 

To first order in T, Eq. (55) gives the generating function as 

w^[j,t 1 ] treB = ij^dt J d 3 x|i^ 0) (x,t)r'(^ 0) (x,t),t) -r(^ 0) (x,t),t) 
- i^ 0) (x,t)r'(^ 0) (x,t),t) +r(^ 0) (x,t),t)|a 3 (t) 

+ ^Jd 3 x<p% ) (x,t 1 )J(x) (59) 

with superscripts (0) and (1) indicating terms of zeroth and first order in Y. 
Eqs. (39) and (40) give the zeroth order fields for general time as 

<^ 0) (x,t) =i J d 3 x' G (x-x',Mi) ^( x . (60) 

( x , t) = -if d 3 x' G$(x - x', t, h) J(x') , (61) 
while from Eq. (41) we find the first-order fields at time t±: 

= JtPx'p dt' |G (x-x / ,t 1 ,t / )a 3 (t / )r / (^ 0) (x / ,t / ),i / ) 

+ G^(x-x',t 1 ,Oa 3 (Or / (^° ) (x',0,t , )| • (62) 

From the definition (36), we see that Go(x — x', t, t') = — Gq(x — x', t', t), 
Eqs. (60)-(62) give the final term in Eq. (59) as 

^ J d 3 x <p%\x,h) J(x) = -^Jd 3 xJ d 3 x' p dt' J(x) 
x |gS(x - x', t>, tt) a 3 (t') F(^ 0) (x', t'),t f ) 
+ G (x - x', t', h) a 3 {t') T'(yf (x', 0, t') } 



so 
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= - \ J rfV dt> I v f (x', h )r' (x', t') , f ) 

This cancels the T' terms in the first two lines of Eq. (59), leaving us with 
the simple first-order result 

W^[JMWo = i J^dt J d 3 x|-r(^ 0) (x,t),t)+r(^ 0) (x,t),t)|a 3 (t) . 

(63) 

According to Eqs. (60) and (61), when we take the functional derivative of 
tp^\y,t) or (pfl)(y,t) with respect to J(x) we get a factor iGo(y — x, t,t\) 
or — iGg(y — x, t,t±), respectively, while the definitions (36) and (49) give 

iG (y-x,Mi) = J d 3 k exp (ik-(y-x)) u k (t)u* k (ti) = (0|/(y, t)/(x, *i)|0> , 

(64) 

-iG5(y-x,Mi) = y d 3 /c exp (ik-(x-y)) u fc (ti)<(t) = (0|/(x, ti)/(y, t)|0) , 

(65) 

which are the same as the factors (47) and (48) associated according to the 
usual graphical rules with external lines. Thus, Eq. (63) is just what is given 
by the "in-in" formalism for a diagram with a single vertex, which can be 
L or R type, and any number n > 3 of external lines for which there is a 
term in T with n field factors. For instance, in the case n = 3, the third 
derivative of Eq. (63) gives 

(0,in|( / 9(x,ti)( / 9(x / ,ti)^(x",ti)|0,in) = i f ' dt f d 3 y r'"(0, t) a 3 (t) 

J —oo J 

x |co(y - x, t, h)G {y - x', t, h)G (y - x", t, h) 

+ G* (y - x, t, h)G* (y - x', t, h)G* (y ~ x", t, ti)| , 

which is the same as given by a graph with a single vertex to which are 
attached three external lines. This three-point function has been calculated 
in a different way, by a perturbative solution for the Heisenberg picture 
fields, by Seery, Malik, and Lyth[6]. It is the analog of the three-point 
function calculated in a more realistic model by Maldacena[3]. 
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Second Order 

So far, we have confirmed in first order that we get the right contribu- 
tions from external lines, but to check that we get the right propagators for 
internal lines, we have to go to second order in T. According to Eq. (55), 
the contribution to W[J, ti]t re e of second order in V is 

^ (2) [J,tl]trcc = i £ MA*) J d3x \- l^R(x,t)r'(^R\x,t),t) 



+ 



^g ) (x,t)^ 1) (x,t)r"(^ 0) (x,t),t) 



2 

+ I^(x,t)r'(^°)(x,t),t) - i^(x,%«(x,t)r"(^°)(x,t),t)| 

+ ±JcPx<p%\x,t 1 ) «/(x) • (66) 
The second-order contribution to </?#(x, fi) is given by Eq. (41) as 

pgWi) = J dV|^^|G (x-x / ,t 1 ,t , )a 3 (t , )r ,/ (^ 0) (x , ,O,i / ) <p { i\ 

+ G* ( X - x', t u f ) a 3 (i') V'^f (x', 0, t>) v4V, f )} • 

By again using the relation Go(x — x', i, i') = — Gq(x — x', i', t) and Eqs. (60)- 
(62), now together with Eq. (67), we see that the final term in Eq. (66) is 

i J d 3 x ^ 2) (x,ti) J(x) = -^Jd 3 xJ d 3 x' p dt' J(x) 

x |gS(x - x', f', *0 a 3 (i') r" (x', 0, t') <p<$ (x', i') 

+ G (x-x',t',t 1 )a 3 (t')r"(^ 0) (x',t'),t')^ 1) (x , ,t')} 

= ~J dV^cZ^j^(x^^ 1 )^''(^(x^^'),^>5 ^ 1) (x^^0 
-4 0) (x',t 1 )r // (4 0) (x',O,t>S 1) x',o}a 3 (O. 

This cancels the T" terms in the first three lines of Eq. (66), leaving us with 
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the much simpler relation 



W^[J,h] trcc = i J^dtaHt) jd 3 x^- ^W^'^W),*) 

+ ^(x,t)r'(^(x,t),t)|. (68) 

(A similar cancelation occurs in each order of perturbation theory.) Eqs. (39) 
and (40) give the first-order contributions to ifL and (p R for general times: 

¥ ,W(x,t) = J d z x' J*' dt' jG(x-x',M')a 3 (<')r'(pf(x',t'),/) 



+ G (x-x',M> d (O 



i*(^V>O,0-^(^*V) 



(69) 



^(x,t) = J d 3 x' p dt' ^G{x-x\t,t')a 3 {t')r'(<pf{x.',t'),t') 



- G5(x-x',M> d (f) 
Using Eqs. (69) and (70) in Eq. (68) gives 

^ (2) [J,ti] tree = - f 1 dta 3 (t) I d 3 x f 1 dt'a 3 (t') I d 3 x' 

4 J — oo J J— OO J 

x | -r'(^ 0) (x,t),t)A i?R (x-x',t,t')r / (^° ) (x / ,t'),i' 

+ r'(v4 0) (x, t),t) a rl (* - x', t, t')r'(v4V, t '),t') 
+ V (x, t),t) a lr (* - x', t, t')r' (v4 0) (x', 0, t') 

- r'^f (x, t), t) a ll (x - x', t, t')r'(v4 0) (x', f 

where 



(70) 



A_ll(x — x', i, i') = — iG(x — x', t, t') + iGo(x — x', t, t') 
A RR (x - x', t, t') = iG(x - x', t, t') - iG* (x - x', t, t') , 
A Li? (x - x', i, t') = iG (x - x', t, t') , 
A fiL (x - x', t, t') = -iG* (x - x', t, t') . 



(71) 



(72) 
(73) 
(74) 
(75) 
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Referring back to the definitions (35) and (36) of the Green's functions, 
and recalling the formula (49) for the interaction picture field, we see that 
Eqs. (72)-(75) give 



A LL (x - x', t, = <0|TMx, tMx.',t')}\0) , 
A RR (x-x',M') = (0|T{ ( p(x,i)¥'(x , ,t / )}|0) , 
A LR (x-x',t,t') = (0|^(x,t)(^(x , ,t , )|0) , 
A RL (x-x',M') = <0| ¥ >(x',t>(M)|0} , 



(76) 
(77) 
(78) 
(79) 



in agreement with the rules (50)- (53) for propagators in the "in-in" for- 
malism. The signs of the four terms in brackets in Eq. (71) are just those 
expected from the rule of associating factors +i and — i with R and L ver- 
tices, respectively. 

We have recovered enough of the results of perturbation theory to assure 
us of the validity of the tree theorem for the "in-in" formalism. 

I am grateful for conversations with P. Greene, E. Komatsu, J. Malda- 
cena, and M. Musso. This material is based upon work supported by the 
National Science Foundation under Grant No. PHY-0455649. 
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